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Q ■ Abstract 
■ 

O ; 

fH , In the hadronic sector of relativistic heavy ion physics, the p ^ 2ti reaction is 

Q^l the strongest one, strong enough to equiUbrate the p with the pions throughout the 

Q^' region from chemical freezeout to thermal freezeout when free-particle interactions 

D . (with no medium-dependent effects) are employed. Above the chiral restoration tem- 

perature, only yo's and vr's are present, in that the chirally restored Ai is equivalent 
to the p and the mesons have an SU(4) symmetry, with no dependence on isospin 
^ ■ and negligible dependence on spin. In the same sense the a and vr are "equivalent" 

scalars. Thus the chirally restored p ^ 2tt exhaust the interspecies transitions. We 
evaluate this reaction at Tc and find it to be much larger than below Tc, certainly 
strong enough to equilibrate the chirally restored mesons just above Tc. When emit- 
ted just below Tc the mesons remain in the Tc + e freezeout distribution, at least 
in the chiral limit because of the Harada-Yamawaki "vector manifestation" that 
requires that mesonic coupling constants go to zero (in the chiral limit) as T goes to 
Tc from below. Our estimates in the chiral limit give deviations in some particle ra- 
tios from the standard scenario (of equilibrium in the hadronic sector just below Tc) 
of about double those indicated experimentally. This may be due to the neglect of 
explicit chiral symmetry breaking in our estimates. We also show that the instanton 
molecules present above Tc are the giant multipole vibrations found by Asakawa, 
Hatsuda and Nakahara and of Wetzorke et al. in lattice gauge calculations. Thus, 
the matter formed by RHIC can equivalently be called: chirally restored mesons, 
instanton molecules, or giant collective vibrations. It is a strongly interacting liquid. 
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1 Introduction 



One of the recurring results from relativistic heavy ion coUisions, emphasized 
by Braun-Munzinger and Stachel and collaborators, found first in AGS and 
then at CERN energies and most recently at RHIC [1] has been the high 
degree of chemical equilibration of the hadronic products of the collisions. All 
of them, with exceptions we shall discuss, freeze out in chemical composition 
at the same temperature, essentially the temperature for chiral restoration Tc- 

Many attempts to explain this equilibration in the accepted scenario of the 
quark gluon plasma(QGP): increased equilibration because of the higher num- 
ber of degrees of freedom in the QGP, or elsewhere, in the greater number of 
degrees of the high temperature hadron gas because of dropping masses, were 
ultimately defeated because of the small couphng constants in these phases. In 
the QGP the perturbative coupling constants were simply too small to bring 
about equilibration. Efforts in the hadronic sector were defeated by the RG 
results of Harada and Yamawaki [2] that hadronic coupling constants go to 
zero at Tc in the chiral limit. 

Recently Brown et al.[3] have shown that a second-order phase transition can 
be constructed for RHIC physics (ignoring the small baryon number in central 
collisions) by making the chirally restored tt and a massless just above Tc. In 
this way there is a continuity with the massless tt and cr, in the chiral limit, 
just below Tc- Quark masses ^ , measured in quenched lattice gauge studies [4] 
are ^ 1 GcV both at |Tc and ST^. We shall assume them to have this value ^ 
right down to T^. In BLRS[3] at T^, the color Coulomb interaction brought the 
TT- and cr-masses down from 2 GeV to 1.5 GeV, and then the 4-point instanton 
molecule interaction brought them the rest of the way to zero. 

The four-point interaction ^ acted as a driving force in an RPA, with all bound 
qq states at unperturbed energy 1.5 GcV in a Furry representation (Coulomb 
eigenstates); i.e., quarkonium. Therefore Brown's [5] degenerate schematic 
model, in which eigenf unctions, etc., have simple analytic forms, could be 
used. 

It was pointed out that the resulting collective excitations, the giant multipole 
states, were seen in the Asakawa et al.[6] LGS, also by Wetzorke et al.[7]. 



The "quark mass" we refer to in this paper is the chirally invariant thermal mass, 
not the current quark mass that breaks chiral symmetry. 

^ As noted in BLRS[3] with our dynamic confinement the mass measured by the 
Polyakov loop is never reached. 

^ The 4-point interaction is modelled as a 5-function with constant coefficient; 
therefore, trivially factorizable. 
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An important new development was the argument by Braun-Munzinger et al. 
[8] that hadron muhiphcities in central high-energy nucleus-nucleus collisions 
are established essentially at the phase boundary between chirally broken and 
chirally restored matter. This must result from multiparticle collisions, which 
are particularly strong just above Tc. 

In this note we develop the schematic model to include nonlinear couplings of 
the giant vibrations just above T^. We study, in particular, the p — > 27r reaction 
in the chirally restored region, both the p and tt being giant resonances. This 
way of looking at the p and tt was developed long ago as the loop sum in the 
Nambu- Jona-Lasinio approach in the chirally broken region of hadrons, but it 
is much simpler here because of the degeneracy in Coulomb-bound qq states 
above T^.. We note that the chirally restored tt's and cr's are the same entities, 
also the chirally restored vectors and axial vectors. Therefore, we need consider 
only the interaction p ^ 27r to determine the nonlinearity of the vibrations, 
this covering all possibilities. 



2 The Free p and tt as Giant Resonances; The Strong p ^ 2tt Tran- 
sition as a Nonlinearity in the Vibrations 

Many research workers are not used to thinking of the p-meson as a giant 
dipole resonance, but that will be a convenient place to begin. The point is 
that the p, all by itself, is a many-body problem in the Nambu- Jona-Lasinio 
language. It is a giant dipole vibration of the negative energy sea, a sum of 
quark- antiquark pieces. 

In NJL at zero density and temperature the negative energy states are filled 
up to a cutoff A ~ 700 MeV. In the modern way of looking at matters, this 
A is the Wilsonian matching scale for constituent quarks [9]. The p or n are 
then obtained as a sum of loop diagrams, the vertex function T equal to 
or 75 for p and tt, respectively, as in random phase approximation, as shown 
in Fig. 1, except that in the chirally broken sector the Coulomb interaction is 
unimportant. The backward going lines represent holes in the negative energy 
sea, the forward-going ones, particles. 

A good complete review of these matters is given by Vogel & Weise [10]. 

Now the width for the free p to decay into two pions is large, ~ 150 MeV, 
meaning that the lifetime of the p is r ~ 4/3 fm/c, still long enough that the p 
can be distinguished as a real particle. The point we wish to make is that the 
p — s> 27r transition in medium can be considered to be a nonlinearity, a giant 
vector vibration changing into two giant pseudoscalar ones. This nonlinear 
coupling will damp the p vibrations, mixing the energy between p and tt ones. 
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Fig. 1. The sum of bubbles in NJL, which can go either forward or backward in 
time, the latter representing ground state correlations. The solid dots represent the 
vertex function, 7^ for the p, 75 for the tt, multiplied by a 5-function and coupling 

constant. The double lines indicate that quark and antiquark, or quark and quark 
holes, are Coulomb bound states; i.e., are in the Furry representation, which we will 
go over to above Tc- 

essentially powering the many-body interactions which give the collectivity to 
these vibrations. 

Note that the p 271 transition is already a strong one in the broken sym- 
metry sector. Also the Ai p + n has an even larger width, and can be 
considered the same way as a collective excitation. However, the enhancement 
due to collective effects, although large, is not nearly as large as it would be if 
the unperturbed energy of all of the particle- antiparticle loops were the same, 
as we shall find to be the case in the chirally restored sector. 



3 Coulomb Bound States Above Tci The Furry Representation 

Shuryak & Zahed [11] introduced Coulomb bound states of quark and an- 
tiquark, or of quark-particle and quark-hole for the region of initial RHIC 
temperatures, suggesting that the breakup of these molecules would help with 
the observed early equihbration. Below T ~ 2Tc the qq states would be bound, 
quite strongly as T moves down towards Tc. The color Coulomb interactions 
can be summed to all orders by going over to quarkonium which has quark 
and antiquark in Coulomb bound states. This is the so called Furry represen- 
tation, well known in atomic physics, and shown by the double lines in Fig. 1. 
This gives us a convenient starting point. 

The cutoff Anjl ~ 700 MeV may be considered as a sort of order-parameter 
for the constituent quarks, as well as Wilsonian matching scale, in that the 
negative energy quarks comprise the condensate which breaks chiral symmetry 
in giving hadron masses. The chiral symmetry breaking scale for hadrons is 
-^xSB — V^Anjl ~ 1 GeV [2]. As T moves above Tc the (collective) chiral 
symmetry breaking state disappears. The chirally restored zero-mass vr, a and 
p become the thermodynamic variables just above T^. This means that the 
scale drops from A^-sb ~ 47r/7r ~ 1 GeV towards the infrared, the scale just 
above Tc being zero. The large distance color Coulomb couphng constant just 
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Fig. 2. The large distance behavior of as{T) from evolution of the Polyakov loop in 
quenched LGS[12]. 

above Tc has been evaluated by F. Zantow [12] from the Polyakov loops in 
quenched LGS. We show the behavior of a{T) in Fig. 2 for T Tc. 

The Bielefeld as{T) does not have the Casimir 4/3 in it that we use. Further- 
more the as is calculated with heavy quarks. In the case of light quarks at 
addition of the velocity-velocity interaction in the appropriate helicity states 
doubles the effective Coulomb interaction. Thus our effective large-distance 
as{T) would reach a value of ~ 16/3 at giving a g^ff ~ 8 for the color 
Coulomb interaction. This is a strong indication of the QCD scale having 
moved far towards the infrared, indicating very low mass modes as thermody- 
namic variables. 

This large distance behavior does not give us the effective color Coulomb 
interaction to be used in calculating our tt, a, p bound states, which are small 
in extent. For these we need the ^^(Tc, r) for small r. This is quite complicated 
because we know that as{T, r) goes to zero as r ^ 0. 

We received from Olaf Kaczmarek [13] results in full (unquenched) QCD. The 
results for the free energy, color singlet internal energy, etc. are very nearly 
the same as in quenched QCD, except that the Tc (quenched) is rescaled to 
Tc (unquenched). We shall assume this to be generally true in what follows, 
although the vibrations have not yet been calculated in full QCD. Note that 
this is full QCD for heavy quarks, so we must add the additional effects for light 
quarks, such as the velocity-velocity interaction and the instanton molecule 
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interaction, which would not enter into the heavy-quark sector. 



We can, however, as in BLRS[3], use the parameter ag obtained from the 
binding of charmonium above T^, namely ccg ~ 0.5. The charmonium atoms 
have about the same radius r ^ 1/3 fm as our a. n. p at T^.. With ag = 0.5, 
doubled to take into account the velocity-velocity interaction together with 
the Nambu-Jona-Lasinio interaction extended up into the chirally restored 
region as done in BGLR[14], BLRS obtained zero energy for their tt and a. 
Below we discuss how the p mass is brought to zero at Tc in the chiral limit. 

We remark that very strong coupling is clearly needed at T = Tc + e because 
the ~ 2 GeV unperturbed energy of the quark and antiquark must be brought 
to zero by the binding potentials in the chiral limit. 

As discussed by BLRS [3], the instanton molecule interactions are expected to 
be unimportant by a temperature T ~ 2Tc. Note that in the chirally restored 
region, the scalar and pseudoscalar degrees of freedom should be the same, 
as well as those of the vector and axial vector. Thus, one can include all 
excitations within the framework of the tt and the p, as we shall do. Further 
note that in accord with Asakawa et al., BLRS found the spin dependence to 
be negligible. 



4 The p ^ 27r Reaction as a Coupling Between Coulomb Bound 



Ultimately we hope to be able to construct the coupling between giant reso- 
nances in the chirally restored region, but as a first step will include only the 
color Coulomb interaction, which leaves out the many-particle aspect of the 
giant resonances. As noted, in the chirally restored region we need consider 
only p's and tt's, which exhaust the chirally restored degrees of freedom. 

In Fig. 3 we show how a qq bound state with the quantum numbers of the p 
can decay into two pions. 

For the transition matrix element (the matrix element to be put into Fermi's 
golden rule) we have 



States. 



M = (^75?(1)^75?(2)5//23?7m?(3)) 



(1) 



where the 6H23 is the Coulomb interaction 




e^rs cos 9 



(2) 
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Fig. 3. How the color flows to show that we make colorless molecules by color 
Coulomb exchange. We show only the single color Coulomb interaction which creates 
an additional quark-antiquark or quark-quark-hole pair. The quark and antiquark 
lines are all Coulomb states. A factor 2 arises because the color Coulomb interaction 
could equally well start from the down-going line on the left; another factor of 2 
because the two tt's could equally contribute as two cr's. 

Here r is the distance between the pion at point 2, which is so small that 
we take it to be a point, and the distributed mass of the p-meson with r.m.s. 
radius [(r|)]^/^ ~ 0.5 fm, at point 3. The distance r will be ~ 1 fm, the average 

distance between particles. We must carry the expansion to obtain the cos^ 
dependence, in order to match to the p which we assume to be polarized in 
the 2;-direction. ^ 



It is of interest to compare the SH of eq. (2) with that of the chirally bro- 
ken p — s> 27r discussed in Sec. 2. We can obtain this from Rcf.[15]. This is 
particularly simple for massless pions ^ which should not be too bad an ap- 
proximation; 

M = g^rup. (3) 
The ratio of the relevant part of eq. (2) to eq. (3) is 



where we have replaced cos by ^1/3, since the square averaged over VL will 

^ In fact, at Tc the p polarized in the time direction is most important. It is related 
to the one in the z-direction by the consistency condition dpfxjdXjx = 0. Thus, one 
can say that they are the same degree of freedom. 
^ We thank M. Prakash for this observation. 
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give 1/3. 



Wc should multiply this ratio by 4 because the Coulomb interaction could 
equally well begin from the antiquark line of the p (left-hand side in Fig. 3) 
and in the chirally restored regime the p can go into two cr's as well as two 
tt's, the TT and a being the same in the chirally restored sector. 

Now, as noted in the Introduction, the radius of the p is ~ 0.5 fm, and we 
found the effective CKg to be ~ 1. We take ~ 5. Since the particles are ~ 1 
fm apart we take r — 1 fm. Note that we do not have asymptotic states in 
the many body calculation, so energy conservation will be somewhat blurred. 
We find R ~ 0.18, i.e. the nonlinearity in the qq Coulomb bound states is 
~ 2/3 as large as that of the chirally broken p when multiplied by the 4 noted 
above. However, the latter is a collective vibration, as calculated in NJL, so 
before we can make a true comparison, we must carry out the same type of 
loop expansion in the chirally restored region, putting in collective effects. 

At Tc, the effective ag is about unity. As will be described in the next section, 
adding loops to the qq bound states and summing them will increase the 
process. Fig. 6 by a factor ~ 4^ = 64 at T = Tf, + e; i.e., the interaction of the 
vibrations at T ^ Tc is much greater than the interaction between qq Coulomb 
bound states. Of course, what we have done is to add up an infinite number 
of attractive contributions in the bubble sum. This is correct in random phase 
approximation when governed by a symmetry principle; i.e., see § 5 Ch. V of 
Brown [5] where the spurious translational state in nuclei is brought to zero 
energy by just such a procedure. The same occurs in the Anderson mode in 
superconductors or the Goldstone mode — our 7r-meson — in particle physics, 
except that these are real, not spurious modes. 

However, in the part of the nuclear many-body problem not governed by sym- 
metry principles, the original Kuo-Brown procedure of including only one bub- 
ble has turned out to be quite good [16]. Of course our situation here is different 
from the nuclear many-body problem, in that we have both the Coulomb in- 
teraction and bubbles. However, the general idea may be the same. In Fig. 4 
we show a typical higher order graph summed by Kirson in his "bubbles in 
bubbles" translated to our Coulomb-|-instanton molecule interactions. 

The "bubbles in bubbles" increase the contributions of repulsive terms - the 
process shown in Fig. 4 is an exchange term - in higher order. Once the 
interactions are more accurately determined; e.g. by LGS, it will be interesting 
to attempt to sum them. (In the case of the nuclear many-body problem, even 
though it is a case of strong interaction, such a systematic treatment has been 
carried out [16].) However, for the moment we will introduce just one bubble, 
assuming the higher-order bubble terms to be cut down by the exchange terms 
with bubbles in bubbles. 
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Fig. 4. Higher-order graph. In this figure we have not used Coulomb bound states 
as representation but bare quark and quark-hole states, in order to make clear the 
connection with the nuclear many-body problem. The dashed lines are the Coulomb 
interaction. We can also attach the bubbles to the quark or quark-hole or to each 
other by instanton molecule interactions. 




Fig. 5. Adding a bubble to one of the tt's in Fig. 3, which also involves adding 
a propagator. The solid dot is the instanton induced interactions. All lines are 
Coulomb bound states. 

Prom the development Fig. 5 and following in BLRS[3] we note that adding a 
bubble as in Fig. 5 will increase the contribution by the factor 



]_GeVF_4 
^-^+1.5GeV"3 

where F ~ 0.5 is the correction for nonlocality. (See Eq. (42) and Appendix 
of BLRS[3]). Now the bubble can be added to either tt or the initial p in Eq. 3 
so that the overall increase is by a factor of (1 -|- 1/3)^ ~ 2.4. 

Including the factor 2.4 in R we find the Golden rule matrix element M to be 
1.6 times that for the chirally broken p ^ 27r, or the width to be a factor of 
(1.6)^ larger, i.e. 
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r = 380 MeV. 



(6) 



Whereas the above is hardly a quantitative calculation, we do end up with a 
width r somewhere in between the nip of 560 MeV of BLRS and the 280 MeV 
we shall end up with here. This is satisfying in that in strong coupling the 
width of a particle should end up more or less equal to its mass. 

In order to come down to the one-bubble insertion as taking care of the role of 
vibrations in the nuclear many-body problem, a lot of calculations, ending in 
Kirson's work[16] had to be carried out. In our present situation wc have both 
color Coulomb, which scales rapidly with energy, and the instanton molecule 
interaction which is tied to the hard glue. Thus a quantitative many-body 
calculation would be difficult. We do, however, have the lattice calculations. 
Dealing with heavy quarks these do not include the velocity-velocity and in- 
stanton molecule interactions, but this may be an advantage for investigating 
how the widths of the resonances grow as T moves downwards towards Tc, 
in that they should not get out of hand, which they do with our estimates 
including everything. 

Given the fact that the coupling must be so strong as to bring the unperturbed 
(chirally invariant) mass 2mq down to zero in order to construct the vr and a 
in the chiral limit, we believe that they will give a strong p 27r transition. 
Since we have only p's and tt's in the chirally restored region above Tc, this will 
result in complete equilibration. (It should be noted that below - but away 
from - Tc the p ^27r reaction is the strongest, and results in equilibration of 
the p with the two pions down to thermal freezeout.) 

We finish this section by showing that the same sort of strong interactions 
that produced the p 27r nonlinearities will give large meson scattering cross 
sections. We draw the color flow in Fig. 6 for pion-pion scattering. In the case 
of the pions, we would use the first term e^/r of 6H, cq. (2) and it is easy to 
see that we will get cross section of hadronic size, since ~ 1 at Tc. 



5 Constructing Vibrations by connecting Coulomb Bound States 

Now the instanton molecule interaction for T ^ Tc is sizable, and effects in 
producing giant collective modes are maximized by the degeneracy in un- 
perturbed energies at the common 1.5 GeV energy of the Coulomb bound 
states. Because of the degeneracy in unperturbed state energy, the problem 
for T > Tc actually resembles more that of the giant dipole state in nuclei 
[17] than NJL for T < Tc. From Ref. [5] one can see that if one includes only 
loops going forward in time, then the energy of the giant collective mode is 
brought down with the entire trace of the secular matrix for these states. The 
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Fig. 6. Pion-pion scattering in the Furry representation. The red-blue and green-red 
contributions to the color Coulomb interaction are shown; of course, the color must 
be summed over. The other color Coulomb (singlet) interactions which establish the 
Furry representation are not shown. 

argument is repeated in Appendix A. Furthermore, if the forward going loops 
bring the energy of the collective mode halfway from the 1.5 GeV unperturbed 
energy to 0.75 GeV, then including the backward going loops ("ground-state 
correlations") will bring the energy all of their way to zero[5]. In other words, 
all of the interactions conspire so as to move the collective mode in energy, 
leaving the remaining non-collective modes at their unperturbed energy of 1.5 
GeV. (As noted earlier, the fact that the pion is a Goldstone mode governs 
the attraction from all sets of bubbles.) 

Indeed, with such strong interactions we will certainly have thermal equilib- 
rium in the region just above Tc (and probably higher). 

Although not yet confirmed by LGS, ^ the scenario by Harada & Yamawaki [2] 
in which the width of the p-meson goes to zero (in the chiral hmit) as T —> Tc 
from below, should be helpful in providing equilibrated hadron emission. In 
their scheme (yfy ~^ as T — > Tc from below. In fact, Brown and Rho [9] 
showed that Qy'^/m*'^ goes smoothly through Tc from LGS of quark number 
susceptibility. Thus if m* 0, gy does also. Furthermore, the G in BGLR[14] 
is essentially the ratio glQq/ml and the rough constancy in G means that as 
goes to zero at Tc in the chiral limit, so does gaQQ- 

Thus, the system of highly equilibrated chirally restored mesons just above 
Tc cut loose at Tc into an environment in which the interactions are zero in 
the chiral limit just below Tc. This is obviously the optimum condition for 
chemically equilibrated particle yields at freezeout, which takes place at Tc, 
as foreseen by Braun-Munzinger et al.[8]. 

^ The necessary unquenched calculations have not been carried out. 
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6 Construction of F — E — to. 



The Helmholtz free energy F — E — ra enters into the exponent of the Wilson 
hne as exp(/3F) and is proportional to the string tension so F = at 
for heavy mesons. In BLRS we found it extremely useful to use the Bielefeld 
lattice results for heavy mesons in charmonium, and to add the additional 
effects relevant for light quarks such as those from the Ampere's law velocity- 
velocity interaction and the instanton molecule interaction. Recently there has 
been a lot of work at Bielefeld evaluating the Helmholtz free energy for heavy 
quarks [18]. It is of interest to calculate F for our chirally restored mesons, 
which add some new aspects, here. 

In the appendix we give the solutions for the instanton molecule Lagrangian 
of BLRS. For any given spin and isospin that the Coulomb bound states are 
coupled to, they form an unperturbed representation degenerate in energy, 
each qq bound state lying at 1.5 GeV. All matrix elements, both diagonal and 
off-diagonal, between unperturbed states are equal, given by the 4-point (zero 
range) instanton molecule interaction constructed by BLRS. For each channel, 
i.e., given spin and isospin, one collective state moves down to zero energy and 
all the other eigenstatcs remain at their unperturbed position of 1.5 GeV. Now 
the partition function is the sum of Boltzmann factors Z{t) = J2n^'^v{~^/'T)- 
For the moment we consider only momentum zero states. The zero-energy 
collective state will contribute unity to the sum. The other states will each 
contribute ~ exp(-1.5GeV/0.175GeV) at T = 7; = 175 MeV, or ~ 2 x 10"^. 

We thus approximate Z for a boson at rest, by 

Z = ^e-^-/^ = l, (7) 

n 

keeping only the zero energy collective state. Thus, the free energy is 

F = -rinZ = 0. (8) 

All that this shows is that in the construction of our chirally restored mesons, 
the large necessary binding energy reduces the large ^ 1 GeV quark masses 
to zero, affecting neither pressure nor entropy, given by 



OF OF 

At temperatures well above T^, where the vibrations are linear, and the epoxy 
is melted - roughly T — 2Tc as will be discussed in the Appendix, we have, to 
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a good approximation, a Boltzmann gas, with each boson carrying e ~ 2. IT 
and (7 = 3.6. The free energy is obviously a minimum. 

However, the connection between equilibration and minimization of the free 
energy is not clear at T^, and possibly not at 1.4Tc, with the vicinity of which 
we are concerned. The usual derivation involves differentiating F to find the 
extremum 

dF^dE- Tda (10) 

and showing that the right hand side is zero in thermal equilibrium. However, 
at Tc the bag constant B, which describes the energy in the epoxy, about 
half of the total bag constant (the soft glue being melted by Tc) increases the 
energy density without changing the entropy. It thus contributes negatively 
to the pressure. For many years lattice gauge simulations tended to find the 
pressure to be negative at which was, of course, unacceptable because the 
system would collapse. Brown ct al.[19] in a crude calculation got the chiral 
restoration temperature to be = 172 McV just by requiring the restored 
phase to begin as soon as the pressure could be made positive. 

Although we don't know whether the Boltzmann gas apphes to the region 

around where the vibrations are very nonlinear, we use it in Sec. 9 to 
discuss the density of the thermal excitations which must be added to bring 
the pressure to zero. 

In the Harada and Yamawaki scheme [2] it is easy to sec why the pressure in the 
hadron gas is low at T^, because the interactions go to zero in the chiral limit. 
(Of course pressures must be equal in both phases in a phase transition.) In 
the resonance gas [20] the pressure on the hadron side is made small by putting 
most of the energy into the binding energy of excited states. 

We can only conclude that the free energy probably is close to zero in the 
chirally restored phase at T^, the value that results for heavy quarks from the 
Wilson loop. 



7 The p-meson is special 

As noted above, because of chiral restoration, we have only two types of chi- 
rally restored mesons tt and p above Tc- The chirally restored Ai is equivalent 
to the p and the chirally restored a is equivalent to the vr. Furthermore, be- 
cause the r in the instanton molecule interaction is a four-vector, and all but 
the last term in the instanton molecule Lagrangian in BLRS[3], to which we 
return in discussion below, involve the square of r, the different isospin states 



13 



p 




p 



% 



Fig. 7. Second order self energy contribution to the p mediated by 2tt exchange. 

are degenerate. In BLRS[3] we showed this to be nearly true for spin effects, 
which were neghgible. Thus, we are left with only the p and n, with their 
various spins and isospins, in an SU(4) multiplet. 

In the chiral limit the vr is constrained as Goldstone boson to have zero mass 
at Tc (just below and just above). Therefore, only the mass of the p can be 
abjudicated; i.e., the only model dependence on meson masses at T^. comes in 
the p-mass. By Tc here we mean coming down to from above; we know that 
going up to Tc below there is a fixed point in the p-mass of zero at Tc [2] . 

Although the mass of the p was found in classical approximation to be zero 
in BLRS[3], the instanton molecule fluctuated quantum mechanically about 
the time axis, with r.m.s. fluctuation in $4, the angle with this axis, of about 

30°, so that rrip ended up at 560 MeV. However, our strong p ^ 271 coupling 
of Sec. 4 means that the lowest-lying eigenmode in the p-channel will be a 
"rhosobar" ; in this coherent linear combination of p and 27r, essentially 



Consider the second-order self energy of the p from this coupling shown in 
Fig. 7 below. Of course this mixes the various 27r states only perturbatively 
with the p which is inadequate for a quantitative estimate. 

If nip is greater than 2m7r, then this loop will involve a principal value integral 
with both positive and negative contributions. Thus, the greatest attraction 
is obtained when rUp = 2771,^. Thus, in the chiral limit rup ~ 0, but with chiral 
symmetry breaking, rUp ~ 280 MeV, twice the pion mass. We do not consider 
this a rigorous argument, but show below that the result is reasonable in 
the sense that going just below Tc where it is emitted the p gets its on-shell 
mass back. This procedure also minimizes the free energy (under the condition 
that the p has the chirally broken mass of 2171,^) which should be achieved for 
equilibrium. We illustrate by this argument that chiral symmetry breaking 
will be very important in determining the mass of the p. Of course this sharp 
nip is used only in lowest-order estimates. With the 380 MeV width we had 
in eq. (6) the p strength will be spread widely up to nearly 500 MeV. 

Note that there is an equal degeneracy of longitudinal p's as pions at Tc. 
Therefore, if the p has a mass of ~ 280 MeV and we neglect the ~ 140 



I rhosobar) = a\p) + Vl — a'^\27r). 



(11) 



14 



MeV difference of tliis mass from tlie pion mass, we woTild ultimately, as tfie 
temperature decreased towards tliermal freezeout, get about as many pions 
from tlie decay of tlie longitudinal p's as from the pions themselves. This would 
be the situation if all interactions cut ofT below Tc, as envisioned by Harada 
and Yamawaki[2] where the width of the p goes to zero in the chiral limit. But, 
of course, this is not so good an approximation, especially not at Tc, where 
our estimates for T^ + e are mixed with those of BSW[8] for Tc — e through the 
explicit breaking of chiral invariance, which converts a phase transition into a 
smooth crossover transition. 

In fact, Braun-Munziugcr et al.[21] note a substantial discrepancy in the po/'K~ 
ratio in semi-central Au-Au collisions. We give their discussion in what follows: 
These mesons have been reconstructed in STAR [22] via their decay channel 
into 2 charged pions. Comparing the preliminary results from STAR with the 
thermal model predictions of BRS[21] , reveals that the measured values exceed 
the calculated values by about a factor of 2. This is unexpected, especially 
considering that BRS use a chemical freezeout temperature of 177 MeV for 
the calculation. One might expect these wide resonances to be formed near to 
thermal freezeout, i.e. at a temperature of about 120 MeV. At this temperature 
the equilibrium value for the Pq/t:" ratio is much smaller than the 0.11 found 
at 177 MeV. Even with a chemical potential for pions of close to the pion 
mass and taking into account the apparent (downwards) mass shift for the 
po it seems difficult to explain the experimentally observed value of about 
0.2. In fact, recently the STAR paper was published [23] where they give the 
final PqI'k' ratio as 0.169 ± 0.003(stat) ± 0.037(syst) for peripheral Au+Au 
collisions, slightly lower than the preliminary ~ 0.2. Other statistical models 
[24,25] also have difficulty in explaining this large a ratio. We return to this 
below. 

Now this excess is difficult to explain with the p width that of the free p, 
V p ~ 150 MeV, during the drop in temperature from 175 to 120 MeV. This 

gives the strongest hadronic interaction and would surely be large enough 
to equilibrate the p's and tt's down to thermal freezeout, the p's changing 
back and forth between two tt's several times. Since they are obviously not 
equilibrated, the assumption that the p-meson has its free-particle width must 
be wrong, and we adduce the Harada and Yamawaki argument to say that 
only that part of the F which originates from explicit chiral symmetry breaking 
remains at Tc\ in any case, that the in- medium width is substantially smaller 
than the free-particle one. 

Let us try the extreme version of the theory in which the p and tt mesons stop 
interacting as they drop below Tc in temperature. In Appendix C we develop 
a schematic model, which has the important ingredients, for how the p of 
mass 2m^ at Tg goes on shell before leaving the system. Out of the po's, only 
the longitudinal one will be of low energy, because of the instanton-molecule 
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polarization, as explained in BLRS[3]. Its strong attractive principal value 
interaction with the pion is maximized in magnitude by a mass of mp = 2771^^ 
(so that there are no negative contributions to the principal value integral). We 
believe the nonzero longitudinal p mass to be the main effect of explicit chiral 
symmetry breaking. We use Boltzmann factors in our schematic estimates. 

The total energy, rest mass plus thermal, of the 770 MeV p at Tc = 175 MeV 
is 1090 MeV, so the Boltzmann factor multiplied by 3 for the spins is 0.006, in 
the standard scenario. In our case only the longitudinal p has the low mass of 
280 MeV at T^. Its total energy is 652 MeV, with Boltzmann factor of 0.024. 
Thus, even though the degeneracy is cut down by 1/3, because of its lower 
mass, its abundance would be 4 times greater. Thereafter, in the chiral limit 
with Tp = 0, the p's would just leave. Given that the explicit chiral symmetry 
breaking is present and that the T(. + e region will be overlapped with the 
Tc — € region in a smooth crossover transition (see the next section) we would 
expect the abundance to be substantially decreased, but still present. The 
experimental enhancement by a factor ~ 2, half of our enhancement in the 
chiral limit, looks reasonable. 

As the longitudinal p goes below T^. in temperature, it will gain back ~ 0.33 
GeV binding energy it got from the color Coulomb interaction just above Tc, 
going on shell. ^ The reconstituted spin interaction also adds to the mass, 
~ 0.25 (mp — m^) = 0.16 GeV. The addition of these to the p-mass of 0.28 
GeV at Tc + e puts the p on shell at Tc — e. 

Although our calculation is made only within SU(2), we believe that the spin 
effects will also be substantially weaker above Tc than below in SU(3). Now 
the Kq-K splitting below Tc is only about half of that of the p-n splitting, 
so the enhancement should not be so strong in the Kg abundance as in that 
of the p, but a factor of ~ 2 in the Kq abundance over that determined by 
equilibration at Tc — e would only improve matters [1], although the difference 
between prediction and experiment is only one standard deviation in Ref.[l]. 

We did not find large deviations from the standard scenario in the STAR 
K*^ / K and (j)/ K*^ ratios. We can see that medium effects will decrease both 
numerator and denominator here. 

In our work we have relied upon chiral invariance to construct the tt and 
a so that their masses go smoothly through chiral restoration. We have not 
considered baryons and don't know how to introduce them. We do know that 
the A cannot have got much of its 300 MeV mass relative to the nucleon from 
the perturbative spin-spin interaction as in the MIT bag model. This would 

We take the below to be that of charmonium, ag ~ 0.33. Just above Tc 
where the Coulomb binding is ~ 0.5 GeV (BLRS[3], Table 1), the effective is 
{as)es ^l.'&o, the 0.33 GeV is 2/3 of the 0.5 GeV. 
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be absent about Tc, with the result that numbers of nucleons and A's would 
be equal there. Whereas STAR shows some excess of A's above the standard 
equilibrium scenario at central rapidity, this excess of ~ 1/3 is small compared 
with the exp(300MeV/Tc) that would be obtained in the perturbative scenario. 
More hkely is that the singlet, isosinglet diquark which is strongly bound, 
and which exists in the nucleon, but not in the A, gives most of the energy 
difference between (possibly incipient) A's and nucleons above T^. The binding 
of the diquark, as our Nambu-Jona Lasinio interaction, goes smoothly up 
through Tc since it results from the 't Hooft instanton interaction. Above Tc 
it will also have an attractive color singlet interaction. Thus, in some way the 
nucleon may be formed above Tc by the diquark picking up a quark, but the 
A would have to be formed by the coalescence of three quarks. In any case, 
the difference in energy between the incipient A's and nucleons just above Tc 
cannot be much smaller than just below Tq- 

Of course, ratios of antiparticle to particle masses will be unchanged by 
whether equilibration is just below or just above Tg. 

There are, thus, lots of ways in which the role of the medium dependence in 
equilibrium is hidden. It is no surprise that its role is clearest in the p/n ratio, 
because most of the p mass is generated dynamically in the chirally broken 
sector, whereas the pion, as a Goldstone boson, remains unaffected (except 
for the explicit chiral symmetry breaking). It is just this behavior of vr and 
p which led to Brown- Rho scaling[27], which was first discovered empirically 
in nuclear spectra, because the vr and p contribute to the tensor force with 
opposite signs. Thus, as the p-mass decreased with density, the p contributed 
more strongly and weakened the tensor force [28] . 

We remark briefly that our above discussion does not include the last term 



<^W = -(V^TmTsV')' (12) 

of the instanton molecule Lagrangian of eq. (37) of BLRS[3]. Above Tc with 
chiral restoration this is the same as 



SLiML = (^7m^)'; (13) 

i.e., attractive coupling of the a;-meson. This may be connected with a discon- 
tinuity in the baryon number chemical potential [29]. 
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8 The Braun-Munzinger, Stachel and Wetterich (BSW) ScenEirio 



Our scenario is that equilibration takes place just above Tc in the chirally 
broken sector. We have constructed a strongly interacting colorless liquid in 
that sector. Nonetheless, our scenario has important points in common with 
that of BSW [8], although they consider equilibration to take place in the 
chirally broken hadronic sector just below Tc. 

These authors make 3 main points: 

(1) Dominance of hadronic reactions with a high number of ingoing particles 
can be realized only very close to the phase transition. 

(2) Two-particle processes are too slow to estabhsh and sustain chemical 
equilibrium near the chemical freezeout temperature T^h- 

(3) Multi-particle scattering is indeed fast enough in order to maintain equi- 
librium for T > Tch- 

These authors carry out an extremely illuminating calculation of the rate 
tq for VL production, through the reaction 27r -|- ?)K — > NVl. As we, they 
take perturbative (Boltzmann) thermal model densities. For the Golden rule 
matrix clement they take the measured p + p — Stt cross section. They check 
by evaluating the rate for production ^ in a semi-classical approach, in which 
the standard two-body rate equation is generalized to multi-particle collisions, 
obtaining essentially the same result. 

Thus BSW reach a time of 2.3 fm, sufficiently short to bring the Q's into 
equilibrium. Close to the VL equilibration scales approximately as cx 

Therefore, in a standard hadronic environment with reasonable parameters, 
chemical equilibration can be determined close to Tc, say at Tc — e. 

In our paper above we have worked at Tc-|-e, i.e., in the chirally restored phase, 
from a different scenario, but with the same general results. Although as we 
formulated it, in the chiral limit, our scenario looks very different from that of 
BSW, we note that introducing the explicit chiral symmetry breaking would 
smooth our sharp chiral restoration transition into a smooth crossover one, 
mixing what happens in our scenario at Tc-|-e to Tc — e. Also, even granted the 
Harada-Yamawaki results in the chiral limit [2], their p- meson width will not 
be zero with the explicit chiral symmetry breaking. Furthermore, with their 
T"^^ factor for the VL~ , equilibration would be likely to continue to Tc — e. 



Because of the triple strangeness of the fi, it may be considered the most difficult 
particle to equiUbrate. 
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Fig. 8. Time tq = riQ/rQ needed to bring 17 baryons into chemical equilibrium via 
multi-particle collisions. The arrow points to the Boltzmann pion density at Tc. 

In practical terms, i.e., fitting the data, the BSW work has been extremely 
successful. Nonetheless, there may be exceptions such as we discuss with the 
p-meson, which cannot be explained in their work. 

We note that inclusion of the + e region in their calculations can only help 
them. At Tc + e we have 8 pions, as compared with 3 below the a, p, 
isosinglet pseudoscalar and pions all being equivalent in the chirally restored 
region. The tq in Fig. 8 goes with the inverse 5th power of the pion density, 
so the higher degeneracy just above would be of help. 

Let us turn the previous statement around and say that the BSW work can 
help us. Namely, with the greater number of pions and the stronger interactions 
the material just above is clearly well and truly equilibrated. 

We believe that our scenario has the following merits: 

(1) Chirally restored particles cannot leave the fireball. Therefore the latter 
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must cool until they can come down to Tc- 
(2) However, equilibration at T = + e and the BSW equilibrium at T = 
Tc — e which increases so strongly with T that even if zero in the chiral 
limit, is hkely to be large because of explicit chiral symmetry breaking, 
will be mixed by the smooth crossover transition, over the region e ~ 5 
MeV, the region below Tc in which BSW find equilibration. Thus, this 
BSW "band" will be the region in which thermally equilibrated mesons 
emerge. 

We beheve that the BSW paper and our note are the beginning of an unrav- 
elling of the puzzle of equilibration. 



9 Discussion 

Brown et al.[19] suggested that the mesonic resonance states tt, cr, p, Ai go 
smoothly through Tc with increasing temperature, simply changing from be- 
ing chirally broken to chirally restored. The same general idea is retained in 
the resonance gas picture [20] in which resonances in addition to the above 
hadrons are included so as to make the transition smooth. As noted in Sec. 2, 
the TT, (J, p, Ai can be constructed as vibrations in the Nambu- Jona-Lasinio 
language below Tc so that our present scenario is not new in making them 
vibrations above Tc- However, it is easier to do the thermodynamics above 
Tc, because the qq molecules form a degenerate unperturbed configuration, 
so that all of the strength ends up in the collective state, none in the non- 
collective ones. This degeneracy endows the ^ Tc region with greatly enhanced 
interaction strengths. 

The pressure must be nonnegative at Tc. This is achieved in the resonance gas 
picture by introducing mesonic excited states as well as resonances, as well as 
other hadrons than we consider. 

Indeed, at Tc we will need additional degrees of freedom because in order 
to bring the p and Ai masses to zero we must polarize them in the time 
direction, and this leaves us with only 8 degrees of freedom, which together 
with the TT and a give 16 at Tc- As noted, with this degeneracy we have only 
~ 1 boson/fm^ at T^ which in Boltzmann approximation gives a pressure 
of p = 0.9Tc/fm^ = 158MeV/fm^, a bit higher when calculated in detail for 
bosons. Now in BGLR[14] we found that ~ half of the ~ 500 MeV bag constant 
was melted in the soft glue by T = Tc, ^ half remaining in the epoxy, or hard 
glue. The —250 McV/fm"' contribution from this to the pressure leaves us with 
a deficit of - 100 MeV/fm^. 

A part of the pressure lost from the transverse p's and Ai's being high in mass 
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will be filled in by pressure from the glueballs. The soft glue that is melted as 
T goes up to Tf. will go into glueballs, but these are high in mass below Tc, 
although certainly not higher than the ~ 1.5 — 2.0 GeV masses of glueballs at 
zero density. The Casimir operator is 3 for gluon pairs, compared with 4/3 for 
quarks, and so the color Coulomb should bring the gg bound states down by 
1.25 GeV. Petreczky et al. find thermal gluon masses to be ~ 20% lower than 
thermal quark masses at T = |Tc and T = 3Tc. With our assumed 1 GeV for 
the thermal quark mass, this would bring the gg bound states down to the 
region of masses where, with a multiplicity of 8 they could bring the pressure 
most of the way to zero. 

Some pressure will come from the kaons: K~^,K',K^,K^ and the vector^ 
K*{892) whose masses at the SU(2)xSU(2) - 175 MeV should be some- 
what decreased below the free masses by the melting of the non-strange quark 
condensates and partial melting of the strange quark condensates. 

Hadrons other than our collective states will enter in and there will be some 
contribution from excited states of our chirally restored n, p, a, Ai although 
these states will he at high energies because of the smallness of our molecules. 

We see no problem in obtaining sufficient pressure well above T^, say at 2Tc, 
because the remainder of the bag constant will go as the epoxy is melted 
and our degrees of freedom will double to 32 as the p and Ai are no longer 
polarized. 



10 Conclusion 

We have argued for equilibration of our excitations just above Tc. Every quark 
and antiquark within the relevant rapidity interval participate equally in a set 
of SU(4) vibrations, the energies of which go to zero in the chiral limit at T^. 
At a higher temperature, approximately that reached in RHIC following the 
color glass stage (t ~ | fm/c), those vibrations are found in the lattice gauge 
calculations of Asakawa et al.[6] and Wetzorke et al.[7]. 

The mesons remain equilibrated, at least in the chiral limit, as they are emit- 
ted, since the interactions have been found by Harada and Yamawaki[2] to go 
to zero as T goes up to T^, from below. And they cannot be emitted until T 
has come down to Tc from above. In the more realistic chirally broken system, 
we believe that chemically equilibrated mesons will emerge in the ~ 5 MeV 
band below Tc of Braun-Munzinger, Stachel and Wetterich [8]. 

^ In fact the vector masses may be considerably less, resembling to some extent the 
/9-meson. 
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Whereas the zero masses of vr and a are protected by chiral symmetry at Tc, 
both shghtly above and shghtly below, the p-meson mass must go on shell at 
770 MeV as it materializes, even though the abundance of p's is determined 
when its mass is nearly zero, slightly above Tc- In the chiral limit this means 
that the number of p's, which are reconstructed from their 27r decay, just below 
Tc should be greater than predicted by Boltzmann factor using their 770 MeV 
on-shell mass. A substantial excess in this number, which appears anomalous 
in the conventional treatment, should remain even after introduction of bare 
quark masses. 

Interactions which dictate the above picture are strong with g ^ 1, very 
different from those in the perturbative quark-gluon plasma often predicted 
for RHIC energies. 

We beheve that we have established, backed by the LGS of Asakawa et al.[6], 
that the matter formed by RHIC is composed of giant collective vibrations of 
qq pairs and powered by the instanton molecule interaction. In a particle- wave 
duality that existed in the work of Nambu-Jona-Lasinio, these vibrations can 
also be interpreted as mesons. Above Tc, these mesons are chirally restored. In 
the formation of these colorless mesons, color is dynamically confined some- 
what above Tc- 

Normally one associates the large and rapid increase with temperature of the 
entropy s in the region of T^, with the formation of the quark gluon plasma. 
However, Koch and Brown [30] (see especially Fig. 3) showed that hadrons 
going massless provided an increase in entropy that matched that found in 
LGS. In fact, somewhat above Tc (unquenched) where the p and Ai are no 
longer polarized in the time direction, we have 32 degrees of freedom in the 
instanton molecule scenario. To the extent that these are essentially massless, 
these 32 bosonic degrees of freedom provide the pressure found in the SU(2) 
lattice gauge calculations [31]. Thus, one might say that the large increase in 
number of degrees of freedom results from the hadron masses going to zero in 
the chiral limit; i.e.. Brown- Rho scaling [27]. 

Thus, RHIC has found a new kind of matter, one we find very interesting. 

In the sense that mesons are composed of quarks and gluons, it may be 
simple-mindedly called "quark-gluon plasma," but hardly so in the sense of 
the predicted perturbative QGP. 

This same new form of matter must have been gone through in the early 
universe, as T decreased from ~ 2Tc to T^. The nice, smooth, essentially 

second-order transition we construct would mitigate against inhomogeneities 
originating in the chiral restoration transition. Since we have dynamical con- 
finement just above Tc, we have not had to discuss deconfinement, which 
presumably requires a discussion of Polyakov lines. 
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A Appendix A: Compcirison with Lattice Gauge Simulation 

As noted in BLRS[3] the main idea of Shuryak and Zahed [11] was that "after 
deconfinement and chiral symmetry restoration at T^., nothing prevents the 
QCD coupling from running to larger values at lower momentum scale until it 
is stopped at the screening mass scale." The lattice scale is ~ (0.5 fm)~^ ~ 400 
MeV because for higher scales (shorter distances) the charges are locked into 
the quarks and antiquarks. In order to compare with the lattice results on as 
shown in Fig. 2, we use the perturbative (long distance) scaling 

asjT^n) _ ln(Ax/AQCD) 
as{T < n) ~ ln(400 MeV/AqcD) 

where we have taken the chiral symmetry scale A^ = 1 GeV and Aqcd = 250 
MeV. The actual increase in Fig. 2 from the as ~ 1/3 for isolated charmonium 
to ~ 8/3 at Tc, where we have included the Casimir operator of 4/3, is more 
than double this, but at least of the right general size, leaving little doubt 
that the increase in color Coulomb at Tc results from movement towards the 
infrared with chiral symmetry restoration. We show below that the instanton 
molecule interaction lowers the scale A essentially to zero. BLRS[3] found that 
the heavy quark Coulomb interaction provided only ~ 1/8 of the attractive 
interaction which brought the n and a masses to zero at T^, so the movement 
towards the infrared is much more pronounced than given by the lattice results, 
which do not include the velocity- velocity interaction and instanton molecule 
interactions, neither of which would contribute appreciably to the heavy quark 
situation. 

The above discussion refers to the large distance behavior. The shorter dis- 
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(A.l) 



tance interaction is strongly influenced by the necessity of the interaction going 
to zero as r ^ 0, because of asymptotic freedom. 

Prehminary results on quenched QCD [13] and in full QCD [32] are given in 
[33]. The chief result of full QCD was to rescale by the Tc (unquenched); i.e., re- 
sults at Tc{quenched) ~ 260 MeV were moved down to Tc{unquenched) ~ 175 
MeV. Whereas we believe this to be true for the heavy-quark color Coulomb 
(singlet) interaction, we also believe that for other problems the instanton 
molecule interaction may introduce an additional scale, as we discuss below. 

The lattice calculations show a great dependence upon radial coordinate r, as 
might be expected from asymptotic freedom and confinement. We showed the 
as from the Polyakov loops in Fig. 2, where as is extremely large at T^. 

However, when all is said and done, the singlet potential energy just above Tc 
can be schematized 



V^y {r>h/2mg) (A.2) 

with as = 0.5 to the accuracy at which we can read off the LGS curves. This 
potential was the one used in BLRS, the Coulomb interaction being modified 
for r < h/2mq, which can be thought of as imposing asymptotic freedom. 
Furthermore, this schematic potential should work reasonably well up to T ~ 
lATc, where the screening mass scale can be read off as ~ 0.5 fm, distances 
below this scale being important in the region of Tc ~ lATc. We believe that 
the above schematization is adequate to discuss physical phenomena in this 
range of energies. Prom T ~ 1.4 — 1.5Tc the Debye screening sets in, so that 
the instanton molecules become bigger and by T ~ 1.9 — 2Tc they become 
unbound. We leave this latter region to Shuryak and Zahed[ll]. 

Whereas we must increase the heavy-quark Coulomb interaction for T ^Tc 
by a factor of ~ 8 in order to include the light quark effects, we have only the 
imderpinning of the heavy quark LGS, aside from our general argument that 
■m„ = rria = at T = Tc in the chiral limit. The LGS show, however, that 
the screening distance of ~ 0.5 fm is not reached until T ~ 1.4Tc, so that the 
charges are still locked into the quark and antiquark up to this temperature to 
a good approximation. This was all that was needed to construct the instanton 
molecule interaction in BLRS, so we assume that there is little change in the 
interactions up to this temperature, although Debye screening sets in rapidly 
above. 

In fact our argument that our scenario at Tc should hold well up to 1.4Tc is 
improved if we construct the collective wave functions. In BLRS we calcu- 
lated the instanton molecule wave functions with only the color Coulomb and 
velocity-velocity potential; i.e., with an ag^eS — 1- This gave us an rms radius 
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Fig. A.l. Gluon condensates taken from Miller (Fig. 2 of BGLR[14]). The lines show 
the trace anomaly for SU(3) denoted by the open circles and the heavier ones by 
filled circles. The Tc marked in the figure is that for quenched QCD, whereas we 
deal with unquenched QCD with Tc = 175 MeV. 

for the TT and cr of ~ 1/3 fm. However, the energy of the tt and a were lowered 
only ~ 1/4 of the way from 2 GcV to 1.5 GeV, as seen in Table 1 of BLRS. We 
then calculated that the instanton molecule interaction would further lower 
the energy to zero, but we did not calculate the wave function of the collective 
state. 

This wave function can be expressed in relative momenta q of quark and 
antiquark. If the energy is to be lowered by 2 GeV, then this relative momenta 
must be made up of components that extend to relative momenta of this size; 
i.e., the coherent wave function, solution of the ^g-bubble sum, must be of 
radius ~ h/2'mqC, or ~ 0.1 fm. In other words, the collective wave function 
made up of all quark- antiquark, or quark, quark-hole components must be 
really tiny if the collective state it is to describe, is to be bound by 2mqC^. 
Thus, the electric charge is really locked into a tiny volume, much less at Tc 
than the screening radius. This is why we believe that there is little change in 
the TT and a masses between Tc and 1.4T'c- 

In fact, we see some indication of the constancy in gluon condensate between 
Tc = 175 MeV and 1.4Tc from the LGS of Miller[34] which we reproduce as our 
Fig. A.l. The initial drop in gluon condensate (trace anomaly) up to Tc ~ 175 
MeV is shown in BGLR[14] to be given by the melting of the soft glue, the 
glue connected with chiral symmetry restoration. It is connected quantitatively 
in BGLR with the energy provided by the dynamical (or constituent) quark 
masses going to zero. The last two points on the right show that the "epoxy" , 
the condensate of hard glue connected with instanton molecules, does not 
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change perceptively from Tc = 175 MeV to lATc- This is just the region we 
discussed above in which the electric charges are trapped in the small instanton 
molecules inside of the screening radii. The large NJL forces, which bring the 
molecule masses to zero or nearly to zero in the region to lAT^ clamp the 
hard glue tightly in place. 

We also got a simple estimate of the molecular breakup temperature from our 
argument that quark and antiquark momenta q must be ~ iriqC. There will be 
break up once the thermal energy of q is equal to rUqC^. But relativistically 
q ~ 3r. Thus T^b ~ rriqC^/?,, roughly 1.97; (unquenched) . 

We now repeat here our treatment of the instanton molecules and their con- 
nection with the giant resonances found in the lattice calculation of Asakawa 
et al.[6] since Ref.[5] is out of print. We have estabhshed that our qq (or quark, 
quark-hole) representation is that of the Coulomb states. Every qq state is con- 
nected to every other qq state by the 4-point instanton molecule interaction 
as in Fig. A. 2. 

We construct this interaction either for quark-particle, quark-hole states or for 
quark- antiquark states, as noted in BLRS. Initially for simplicity we consider 
only states of total momentum zero; i.e., the qq state has no translational 
energy. The latter will be put in later. 

As in BLRS we take the quark and antiquark masses to be 1 GeV, considering 
all qq state energies initially at 2 GeV to be lowered to 1.5 GeV by the color 
Coulomb interaction. The diagonalized collective state will turn out to be a 
boson state at zero energy and zero momentum. We will then treat this state 
as a boson, neglecting the incoherent states which remain at 1.5 GeV, where 
their Boltzmann factors at Tc are so small that they can be neglected. 

The instanton molecule interaction starts out as a nonlocal one, the nonlocality 
being governed by the of the instanton zero modes. The is sharply 
peaked, mostly lying within a radius of r ~ ^j2/bp, where p ~ 1/3 fm is the 
radius of the instanton (Appendix of BLRS [3]). This regulates the nonlocality 
and acts as a cutoff in the possible momentum range of the qq states. The way 
this is handled in BLRS is to use a 5-function 4-point interaction, but to cut 
it down by a factor F = (0.75)^. The underlying nonlocahty restricts us to a 
band in momentum space. 

We have, therefore, several approximations, but we are guided in our numerics 
by the principle that the pion energy must turn out to be zero in the chiral 
limit in order to make a smooth chiral restoration transition. Of course, the 
chirally restored cr-meson must also have zero mass at T = Tg. 

If we keep only qq states going forward in time (ignore ground-state correla- 
tions) we have the quantum mechanical problem to solve, for T — Tc, 
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Fig. A. 2. Instanton molecule interaction g which connects qq bound states. 
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(A.3) 



In other words, having reduced the instanton molecule interaction to a 4- 
point one, all matrix elements in the secular matrix are equal. Here gi is the 
instanton molecule interaction. The first —0.5 GeV times the A^-dimensional 
unit matrix Inxn is just the color Coulomb interaction. 

Now from the principle of insufficient reason (all off-diagonal matrix elements 
are equal) we know that one eigenfunction (the collective one) is 



V'coU = 



/ 1 \ 



(A.4) 



V 1 / 

and we easily find that 



N 



E = -0.5 GeV - N gj. 



(A.5) 



In BLRS[3] wc accomplished the same diagonalization by a sum of qq loops 
going forward in time, giving 



Ngi = 0.75 GeV. 



(A.6) 



Thus, at this level the collective state is brought down from the 2mq — 2 GeV 
to 0.75 MeV. 

The equivalence between solving the secular matrix and summing loops as 
was done in BLRS going forward in time was developed in detail by Brown 
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Fig. A. 3. Backward-going loops; equivalently, ground-state correlations. 
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Fig. A. 4. Spectral functions of Asakawa et al.[6]. Left panel: for = 54 (T ~ lAT^). 
Right panel: for iV^ = 40 (T ~ l.OTc). 

[5]. There it was shown that if the collective state is brought down half-way 
to zero (In our above case from 1.5 GeV to 0.75 GeV by forward going loops.) 
then with inclusion of backward going ones it will be brought to E = 0. The 
backward-going diagrams are really ground-state correlations, are shown in 
Fig. A.3. 

We consider only the fundamental mode of the vibrations, the higher modes 
being lattice artifacts [35,36]. 

Let us consider the left panel of Fig. A. 4, taking the fundamental to be at 
2.1 GeV. As noted, BLRS[3] had no isospin dependence and very small spin 
dependence, so the SU(4) character of the vibrations is explainable. Now we 
developed above that the situation at T = 1.4Tc should not have changed much 
since Tc, because the screening scale has only reached that of the instanton 
molecules about there. 

In fact, the calculation in BLRS[3] used the Coulomb wave functions and 
then summed the loops with four-point interaction given by the instanton 
molecule one. The very small collective wave function of size r ~ h/2mq was 
not constructed. But the Coulomb potential was chosen so that dV/dr\r=o = 
and the instanton molecule source ^jJ■ip ~ (r^ -|- p^)^^"^ also had dV/dr = 
0, both of which should follow from asymptotic freedom. Furthermore, since 
the molecule energy is a minimum at T^, it should recover (from zero) only 
quadratically as (T — T^)^ from its minimum. All of these effects make the 
growth in molecule slow as T moves above T^, so that situation at T = 1.4Tc 
should not be much different from T = 0. 
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Neglecting small changes from Tc to lATc in nia and m^r, we find that with 
rriq = 1.2 GeV, giving unperturbed molecule energy of rrig + rrig = 2.4 GeV, 
then the heavy quark Coulomb potential will bring this down ~ 1/8 of the 
way to zero, or to 2.1 GeV. This niq is ~ 20% below the Petreczky et al. 1.5 
GeV mass from quenched calculation at 3/2 but still above the 1 GeV 
that would be applicable if rriq scaled with Tc in going from quenched to 
unquenched calculations. Although this difference is probably not larger than 
uncertainties in the LGS, we believe that it is indicative of another scale than 
the color Goulomb one, that of the instanton molecules, between Tc and lATc. 

In any case, the fact that the LGS find a vibration of energy hcu ^ 2 GeV 
at T = lATc and that the vibration requires attraction between quark and 
antiquark, shows that rUg > 1 GeV. 

Now using the fact that the LGS have only ~ 1/8 of the total attraction, 
we find that the rUfj and ttItt would still be essentially zero at T = lAT^. We 
confess that our above scenario seems to fit a complicated situation too well, 
but believe it to be roughly true. Notice that at T = 1.4Tc, where we have 
argued the situation is not much different from T^. because the charges are 
locked into molecules inside of the Debye screening, the LGS show a complete 
SU(4) symmetry. Harada and Yamawaki[2] show that in the chiral limit the p 
mass joins the tt and a mass as T goes upwards to Tc at zero. It seems most 
reasonable that the p mass is also zero just above Tc. 

Once the screening range of ~ 0.5 fm is reached by the size of the instan- 
ton molecules, the latter would be expected to become rapidly unimportant, 
probably by T = 1.9Tc where the situation appears to have gone perturbative 
[3] or at least to have reached the breakup of the qq bound states. Since the 
instanton molecules are unimportant at 1.9Tc, that temperature should scale 
with Tc in going from quenched to unquenched; thus in full QCD it should be 
1.9Tc{unquenched) ~ 332 GeV, just above the temperature at RHIC following 
the color glass stage. 

In fact we referred to the LGS as having seen vibrations. From our above 
discussion it should be clear that they saw only the qq bound states, and only 
the binding of these by the Coulomb interaction, not including the velocity- 
velocity interaction or instanton molecule interaction which do not contribute 
for heavy quarks. The real vibrations are those we calculated earlier in this 
section, with inclusion of backward going graphs (ground-state correlations). 
In order to get the full collectivity these latter effects have to be included, 
but they cannot be in Euclidean time LGS. However, they can be added 
theoretically as we did for T — Tc. 

Thus, for the modes such as the tt and a which go massless at Tc, the best that 
even an unquenched calculation (full QCD for light quarks) can do is to bring 
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them down halfway to zero because the ground state correlations cannot be 
handled in Euclidean space. (An exponentially decreasing antiquark cannot be 
a quark going backward in time.) (We noted above that in Minkowski space, 
if the forward going graphs brought the vibrations down halfway, then the 
inclusion of backward-going graphs would bring them down all the way. Note 
that our schematic mode eqs. (A.3)-(A.5) is equivalent to a treatment with 
only forward-going graphs.) In this case, because the diagonalization of the 
secular matrix puts all of the trace into the one collective state for each spin 
and isospin, the strength in this one collective state is equal to the sum of 
the strengths in all of the individual ^g-molecules, although the energy of the 
collective state will be moved down 4-times further at Tc. 

The theory of vibrations has many names; e.g., random phase approximation, 
particle-hole solutions of the Bethe-Salpeter equation, bosonization of particle- 
hole excitations, time-dependent Hartree-Fock, linearization of the equation of 
motion. The latter carries the best description; they are the best that can be 
made from linearization. They are also the lowest-order approximation that 
satisfies conservation equations; e.g., conserves Newton's law[37]. 

It may seem somewhat amazing that we can apply the same theory of vi- 
brations in matter at temperatures of ~ 200 MeV as we do in nuclei at low 
energies. However by eliminating longitudinal and scalar gluons in favor of an 
instantaneous Coulomb interaction which gives the Furry interaction, and 
then, integrating over time and regularizing by integrating over the space- 
like nonlocality, one gets an ~ 50% reduced, but instantaneous interaction 
(Appendix of BLRS [3]). Thus, we have the same tools in hand as in low- 
energy nuclear physics with the added simplicity that our Furry representa- 
tion quark-antiquark states are really degenerate, as confirmed by the LGS. 
(In low-energy nuclear physics the term energy in the giant dipole excitation 
huj is the distance between shells but the unperturbed energies are split by 
the spin-orbit interaction which gives added structure [17].) Thus, we do have 
the same problem, but heavy quark LGS have seen only a small part of it, as 
we laid out above. Nonetheless, the small parts of the "iceberg" seen in the 
LGS give already the most prominent excitations above Tc. 

There may, in fact, be an advantage in only a small part of the light-quark 
interaction entering into the LGS. The vibrations look pretty linear at T = 
lATf,. However, with the ~ 8 times larger interaction and the strong p ^ 2^ 
mixing we have investigated, the widths of the p and tt might have been 
spread out over such a wide range of energies that the vibrations would not 
be discernible from the background. 



" Note the movement towards the strong = 4 supersymmetric Yang-Mills theory 
at finite temperature in which the scalar and longitudinal phases travel faster at 
superluminal speeds as the coupling becomes stronger [38]. 
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B Appendix B: Consequences of Rescaling 

We needed Appendix A in order to discuss the physics accompanying rescahng, 
primarily to clarify that the LGS takes note of the change in scale from A^^ = 
1 GeV to the instanton molecule (or, equivalently, chirally restored meson) 
scale set by their zero masses, in the chiral hmit. Thus, ahhough we go up 
in temperature to Tc, suddenly the scale governing the behavior of our 
thermodynamic quantities drops out, and the system reverts way back toward 
the infrared. However, until chiral symmetry is restored, the system knows 
nothing about the infrared behavior that enters at Tc. 

Thus, within the chirally broken system there is a movement towards weak 
coupling in the effective variables, i.e., the hadrons, resembling asymptotic 
freedom. At the risk of a misnomer (since non- asymptotic fixed point is in- 
volved), we shall call this "effective-sector asymptotic freedom." Beta func- 
tions are written down in these variables, the vector mesons replacing gluons, 
etc. in the hidden local symmetry of Harada and Yamawaki[2]. The flows of 
various quantities in their renormalization group treatment join smoothly with 
those of QCD, so their theory in the effective sector mimics QCD, even to the 
point of having a local gauge symmetry: The critical point T^. delineates flavor 
gauge theory from color gauge theory, the latter ceding to the former in a 
smooth way . One may say that the HLS theory is the shadow of QCD in 
the effective sector. It even has its own gluon condensate [14], the soft glue 
which melts as the temperature increases towards Tc (see fig. A.l) although 
the effective sector knows nothing of the hard glue, the explicitly broken chiral 
symmetry above Tc. Thus, in their own effective world the effective theory has 
a gauge theory description, effective asymptotic freedom, etc. It gets its scale 
from the spontaneous breaking of chiral symmetry through the condensate of 
soft glue [27] as shown in Fig. A.l. 

In fact, from the region in which T ~ 125 MeV up to Tc = 175 MeV the 
system behaves as if it is in a mixed phase with nearly constant pressure since 
as described in [14] the increase in energy with temperature goes into melting 
the soft glue. Thus the elliptic fiow is very low and, indeed, collapses around 
mid-rapidity for protons in the NA49 40 AGeV experiments [39]; the velocity 
of sound V = p/E decreases with increasing energy. It is in this region that 
the Harada and Yamawaki[2] scenario of vector manifestation brings out the 
effective sector asymptotic freedom. The hadrons are mainly pions and vector 
mesons and the fixed point of nip — and gv — ets Tc is reached from below 
comes in to further soften the interactions. In BGLR [14] we showed that 
the nucleons dissolved into constituent quarks as the soft glue is melted, the 



This "continuity" of gauge degrees of freedom between color gauge and flavor 
gauge was already conjectured in [40]. 
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quarks then going massless, in the chiral hmit, as T moves up to Tc. As noted 
at several places in our paper, the pressure is approximately, possibly exactly, 
zero at Tc- (The free energy F is essentially —pV for our purpose, so making 
it zero makes the pressure zero.) All of this is efficiently encoded in the HLS 
with the vector manifestation, HLS/VM, i.e., the effective gauge theory in 
effective variables. Whereas it is true that this theory thus far describes only 
the behavior of pions and vector mesons, and has not been extended to include 
baryons , the former are the important variables for the thermodynamics 
of the relativistic heavy ion coUisions. 

Now the effective sector given by the Harada-Yamawaki theory which is not 
valid above knows nothing about the sector in which chiral symmetry is 
restored. Before Harada and Yamawaki RG[2], it was clear that something 
must change, because the Hagedorn temperature limited the increase in Tc- 
What we have shown is that simply to enforce the most basic (and simple) 
requirement that the Goldstone boson, the pion, move up smoothly through 
Tc without changing mass, requires the weakly interacting effective sector at 
Tc — e to be replaced by the most strongly interacting sector in the gamut of 
the system, and this is what RHIC has done. 

The "new" system, representing bona- fide QCD, has been greatly studied and 
described. It has its own gluon condensate, the hard glue or epoxy, which ex- 
plicitly breaks chiral symmetry. QCD has its own running coupling constants, 
so the reincarnation of asymptotic freedom, the shadow of which occurred just 
below Tc in the behavior of the effective variables, occurs only higher up, at 
T ^ 2Tc, possibly above the maximum temperature reached at RHIC. 

We thus see that our rescaling arguments reverse the generally accepted sce- 
nario. Coing up in temperature we first encounter a weakly interacting region 
in the effective variables, just below Tc. With chiral restoration this gives way 
to the most strongly interacting matter encountered in heavy ion collisions. 

With the reversal of the commonly accepted scenario, the weakly interacting 
system below Tc changing into the strongly interacting system just above Tc, 

it is clear that equilibrium has to take place in the latter region. We have 
adduced arguments from LGS and from experiment (especially the p-meson 
decay) to support this new and surprising scenario. 



It has been shown however that the VM holds also for the constituent quarks, so 
the dynamical quark mass vanishes at the chiral transition [41]. 

Recall that the pion mass is brought down from mq + rriq ~ 2 GeV to zero by the 
interactions! 
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C Appendix C: A Schematic Model for How the p Goes On Shell 



Shuryak and Brown [42] calculated that the p-meson freezes out at a temper- 
ature of 120 MeV and a baryon density of 15% nuclear matter density. The 
p mass at this density was measured by STAR [22] to be 700 MeV. The p 
then decoupled from the system and got the remaining ~ 10% of its on-shell 
mass back from its kinetic energy. In the regime below 125 MeV density de- 
pendent effects, which we do not discuss here, dominate, but from 125 MeV 
up to Tc = 175 MeV thermal effects connected with the melting of the soft 
glue predominate [14]. In this paper the constituent quark mass was correlated 
with the binding energy of the soft glue (which is responsible for the dynam- 
ical breaking of scale invariance; i.e. Brown-Rho scaling. The soft glue was 
shown in Fig. A.l. Shuryak and Brown (see Sec. II. E of Ref. [42]) argued that 
collision broadening built the F* back up to the on-shell Fp = 150 MeV at 
Tfreezeout — 120 McV; the part of the width connected with the two-7r decay 
being brought down to 100 MeV because of the reduced (p-wave) penetrability 
involved in the decay. We consider here only the p connection with the two-7r 
system. 

In the Harada and Yamawaki[2] scenario m* and F* go to zero in the chiral 
limit as T goes to T,, from below; also gv/m* goes as a constant as T goes to 
Tp. This means that ^'^ 




F 

three powers coming from the p-wave penetrability and two from {gyY 



(C.l) 



We now make the assumption that the dynamically generated part of the 
p-mass scales with the constituent quark mass in the chiral limit. This is a 
reasonable assumption near as it was near ric [41], but we extrapolate it all 



Note that from Sec. II. E of Shuryak and Brown[42], r*/rp would go to zero as 
(m*/mp)^ even if gy did not scale, because of the p-wave pion penetrabilities and 
reduced phase space for decay. Indeed the expression for the decay of the free p is 

" V47ry 12 V 

If we replace by m* and let the 700 MeV nearly on-shell mass found at thermal 

freezeout by STAR[22] go to zero, this expression gives a F* that scales rapidly 
dropping to 44 McV by the time m* = 350 MeV. Thus, our decoupling of the p 
from the pions does not depend much upon the scaling of gy \ matters would not be 
much changed if one had only BR scaling. 
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Table C.l 

Effective decay widtlis and mean free paths for various temperatures. 

of the way until the constituent quark gets nearly all of its mass back. From 
BGLR[14] we have the binding energy of the soft glue going as 



B.E. = 




We recall that this equation is motivated by Nambu-Jona-Lasinio in which the 
binding energy in the vacuum - in this case that of the soft glue - is obtained 
in terms of the negative energy quarks in the condensate which breaks chiral 
symmetry having dynamically generated masses rriQ. 

For rriQ <^ k we get from (C.2) 



B.E. 
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(C.3) 



We then make the simplest schematic model incorporating the above ideas as 



r* w 100 MeV 

which comes to 100 MeV at T=125 MeV rather than at 120 MeV as found 
by Shuryak and Brown[42], but we want a point midway between 120 and 
130 MeV to do our later integration as the STAR experiment [22] does not 
constrain F accurately. (In fact, there are indications that at densities < Uq 
and probably also at low temperatures - see Fig. 2 of Brown and Rho [43] - 
g* does not scale as rapidly as m*. This would move the point where F* = 100 
MeV lower.) We then find the results in Table C.l, where (F*)"-*^ gives an 
indication of the mean free path in fermis. 

As noted in BGLR[14], this region of temperatures is in the environment usu- 
ally described as a mixed phase. Rather than being at constant temperature 



1 - 



700 MeV y V50 MeV 



+ 



700 MeV. 



(C.4) 
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Tc — 175 MeV, where the bag constant is usually assumed to be a source of en- 
ergy going down in temperature, so that it keeps the temperature constant in 
a real mixed phase, in our case (as shown in the lattice calculations. Fig. A.l) 
the energy is furnished over a relatively large (175 - 125 MeV) range of tem- 
peratures, giving an ~ 50 MeV range of mixed phase. The usually assumed 
~ 5 fm/c from phase transition to freeze out, means that each interval of 10 
MeV takes the system ~ 1 fm/c to traverse. We see that the p mesons are 
never really equilibrated, although their final decay rate just before freezeout 
would be sufficient to equihbrate them had it acted for the full ~ 5 fm/c. They 
mostly decay in the last 2 fm/c, centered about T — 130 MeV; i.e., from 140 
to 120 MeV, the latter the freezeout temperature. In that interval, by taking 
the average mean free path as 2.5 fm, roughly exp(— 2fm/2.5fm) = 0.45 of 
the p mesons don't decay into two pions. Now the pions are always in equi- 
librium (with themselves). Thus, the rate at which pions create p mesons in 
the last 2 fm/c would be equal to the rate at which p mesons go into pions, 
were there equilibrium. But the pion equilibration doesn't depend apprecia- 
bly on the p mesons, so we can use the equilibrium po/7i~ = 0.08 [24]; i.e., 
(1 — exp(— 2/2.5)) X (po/tt") = 0.55 x 0.08 will be created from the pions in 
the last 2 fm/c. (This reverse process uses the same F* as the p-decay.) Thus, 
at thermal freezeout we have 0. 45 (pq/tt") initial left, plus the 0.55 x 0.08(po/7i'~) 
produced in the last 2 fm/c from the pions, the sum of which we equate with 
the STAR result 

0.45 X (po/7r~)i„itiai + 0.55 x 0.08 = 0.169 ± 0.037 (C.5) 
which gives 

(po/7r-)initiai = 0.20 - 0.36. (C.6) 

The systematic error of ±0.037 gives the rather large uncertainty in the initial 
number of p's. We come to the factor of ~ 2 — 4 increase over the Braum- 
Munzinger et al's 0.11 value of (pq/tt^) equilibrated at 177 MeV assumed just 
below Tc. In the text we used m* = 2m^ just above so as to minimize the 
free energy, but with out estimated F* ~ 380 MeV just above Tc the m* will 
be spread over a broad peak lying somewhat below the on-shell p mass of 770 
MeV. 

We believe our schematic model illustrates the main physics: 

(i) The explicit chiral symmetry breaking is important, giving the 2171^^ in 
the p-mass which roughly gives the center of the p mass distribution. 

(ii) The interaction of the p with the pions is negligible at Tc, increasing as 
T drops, and becoming large at the end of the mixed phase as the p is 
nearly on shell. 
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Fig. C.l. Dilepton production through the p. 

(iii) The initial abundance of p's is large, because they are equilibrated just 
above Tc with a spread of masses well below the mass of the on-shell p. 

Given (iii) it is very reasonable that some of this initial abundance remains at 

thermal frcczcout. In fact, we believe in the lower limit of at least twice that 
of the standard scenario in which the p freezes out just below Tg, because with 
large pion and p chemical potentials it seems possible to nearly preserve the 
factor 2 from chemical to thermal freezeout. 

As suggested before, CERN, in finding a weakly interacting system, was oper- 
ating mainly in the mixed phase described in Appendix B. RHIC clearly comes 
down through the region of temperatures well above Tc, through lATc down to 
Tc where we have our low- mass chirally restored mesons. In this ~ 2 fm/c the 
p-mesons have time to decay into dileptons. This should substantially increase 
the invariant masses in our range up to ~ 400 MeV and possibly higher. 

We now go on to discuss how our scenario of equilibrium above Tc will affect 
the dileptons. 

Vector dominance is known to be violated due to the vector manifestation fixed 
point a = 1 at finite density and/or temperature [2]. However the "intrinsic" 
violation more prominent in the baryonic sector [9] is not expected to make a 
qualitative influence in the mesonic sector, so we will simply adopt the vector 
dominance in our discussions. We need however take into account that the 
vector operating on the physical (finite density) vacuum does not only create 
the vector meson (and 7-ray), but there are other many-body vector excita- 
tions, especially the Rapp/Wambach "rhosobar." The latter is an A^*(1520) 
excitation coupled to a nuclcon hole, with vector quantum numbers. In fusing 
Brown/Rho and Rapp/Wambach [43] Rapp considered the production of rho 
mesons by pions, with the later decay of the p into dileptons in Fig. C.l. The 
TTTrp coupling involves a Qy and the p-propagator, (/c^ + m*^)~^. Thus, the 
vector dominance coupling at the photon point A;^ = is given by 



em; 
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(C.7) 
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Fig. C.2. Dilepton production via rhosobar dominance. The rhosobar is subsumed 
in the big dot on the left. 

for gauge invariance. Now as we go away from the photon point to produce 
dileptons, in general the invariant mass of the latter is ~ m*, the k of the 
photon being small compared with the /cq. (Back-to-back kinematics involves 
only an ~ 10% error in dilepton production.) In doing this, we extend the 
vector dominance to the rhosobar, treated as a vector particle. 

Thus, using vector dominance. Fig. C.l reduces to Fig. C.2, where the p being 
subsumed in the big dot on the left. 

We can check the validity of this approach at zero density and temperature 
from the particle-data book where Fp = 150 MeV and the branching ratio into 
e+e- is (4.49 ± 0.22) x 10-^ or 

Fe^ = 6.7keV. (C.8) 

In Fig. C.2 the branching ratio would be simply = 5.3 x 10^^, a difference 
of only ~ 10%. 

We see from Tab. C.l that F* is large only in the last part of the mixed phase, 
from T = 130 to 120 MeV, the nearly hadron phase as the p-meson goes back 
on shell. We will call this the hadronic phase because the hadrons are nearly 
on shell. We assume the Fg^ = a^F*, so that dileptons come only from the 
hadronic phase in CERES. 

Rapp and Wambach[44] enhance the dileptons substantially by introducing 
the rhosobar, which at zero temperature puts ~ 20% of the p strength at 
590 McV. With temperature and many-body effects the width of the rhosobar 
broadens to F ~ 250 MeV, so that strength is moved down another ~ 125 
MeV. With tails of strength functions, a substantial number of dileptons are 
produced down to ~ 250 MeV. 

The factor m*^ / gy must be introduced into the amplitude for dilepton pro- 
duction if Brown- Rho scaling is to be "fused" with Rapp/Wambach. As noted 
by Brown and Rho [43] this factor cancels the enhancement that would be ob- 
tained because of the B-R dropping p-mass, so that Rapp/Wambach should 
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Fig. C.3. We redraw Fig. C.l in the chirally restored region, tt+tt ^ p ^ 7 — e^e , 
in the one-bubble addition we used to obtain the T = 380 MeV of eq. (6). 

fit well the hadronic dileptons. 

We go into some detail with this because we believe that the hadronic dilep- 
tons in RHIC will be similar to those from the 200 GeV/nucleon CERES 
experiments. The scalar density in RHIC is not much less than at CERN be- 
cause the antibaryons give a large contribution in the former. The factor of 
2-4 in the p abundance at Tc will enhance the RHIC dileptons. However, the 
chirally restored region of temperatures will contribute at RHIC. In particu- 
lar, the region of temperatures from to ~ 1.4Tc will give dileptons in the 
range from to ~ 400 MeV in our scenario, the broad region being given by 
the large width [6] . The contribution to the higher dilepton invariant masses 
from higher temperatures will be smaller because of the lower densities at the 
origin for these vibrations, this density relevant for dilepton production. 

The phenomenological vector dominance implies that a strong interaction, in 
Fig. C.2 the p-mcson hidden in the large black dot, precede a photon in the 
electromagnetic decay. Since, as noted in BLRS[3] the vector mesons move 
smoothly up through Tc, we believe that this assumption remains valid above 
Tc for the chirally restored mesons. In fact, the chirally restored p is the only 
vector particle above Tc, the rhosobar having disappeared with the nucleons. 
In any case, it seems the simplest way to couple dileptons to the p; namely, to 
add an off-shell 7-ray to the p (which is itself a sum of bubbles in the BLRS 
random phase description) as in Fig. C.l. This gives us the simple estimate 
for the chirally restored sector of 



Tcm ^ a'r, (C.9) 

where we have estimated Tp to be ~ 380 MeV a,t Tc + e in eq. (6). 

One can see that eq. (C.9) is consistent with our calculation giving eq. (6) of 
r where we have added one bubble. Fig. 5, to each of the pions and p in fig. 3. 
We redraw the final tt's and p in Fig. C.3, to recover the analog of Fig. C.l in 
the chirally restored region. 

We thus expect about equal dilepton abundances from the chirally restored 
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and hardonic sectors, with very little contribution from the mixed phase. The 
hadronic abundance should be about the same as in RHIC, and that from 
the chirally restored sector to be at lower energies, centered about 2m^, but 
with a large width which causes the upper end to overlap with the hadronic 
dileptons. 

Because of the very strong coupling in the chirally restored sector for T ~ 
Tc — lATc, we cannot make any really quantitative calculations, but we believe 
our schematic model to give the main features. It will be exciting to confront 
them with experiment. 

Of course, we know that the dilepton "cocktail" (background) is more than an 
order of magnitude greater in the region about dilepton invariant mass ~ 2m.^, 
so the low mass dileptons from the chirally restored sector may be difficult to 
separate from these, but hopefully there will be distinguishing characteristics. 



References 



[1] P. Braun-Munzinger, D. Magestro, K. Redlich, and J. Stachel, Phys. Lett. B 
518 (2001) 41. 

[2] M. Harada and K. Yamawaki, Phys. Rep. 381 (2003) 1. 

[3] G.E. Brown, C.-H. Lee, M. Rho and E.V. Shuryak, hep-ph/ MmTSj Nucl. Phys. 
A, accepted, (denoted as BLRS) 

[4] P. Petreczky, F. Karsch, E. Laermann, S. Stickan and I. Wetzorke, Nucl. Phys. 
Proc. SuppL 106 (2002) 513; |hep-lat/0110111 



[5] G.E. Brown, "Unified Theory of Nucleon Models and Forces'^ 1967, North 
Holland Pub. Co., Amsterdam. 

[6] M. Asakawa, T. Hatsuda and Y. Nakahara, Nucl. Phys. A715 (2003) 863c. 

[7] I. Wetzorke, F. Karsch, E. Laermann, P. Petreczky and S. Sickan, Nucl. Phys. 
Proc. Suppl. 106 (2002) 510. 

[8] P. Braun-Munzinger, J. Stachel and C. Wetterich, |nucl-th/0311005| 

[9] G.E. Brown and M. Rho, |nucl-th/0206021[ Phys. Repts., in press. 

[10] U. Vogel and W. Weise, "Progress in Particle and Nuclear Physics" 27 (1991) 
195. 



[11] E. Shuryak and I. Zahed, hep-ph/0307267 



[12] F. Zantow, 2004 University of Bielefeld Thesis, unpublished. 

[13] O. Kaczmarek, F. Karsch, P. Petreczky and F. Zantow, ^hep-lat /0309121| 



39 



[14 

[15 

[16; 
[is: 

[19 

[2o; 

[21 

[22; 
[23; 
[24 
[25 
[26; 
[27; 
[28; 
[29; 

[3o; 

[31 

[32; 
[33; 

[34; 

[35; 
[36; 
[37; 



G.E. Brown, L. Grandchamp, C.-H. Lee and M. Rho, Physics Reports, 391 
(2004) 353. (denoted as BGLR). 

"Field Theory in Particle Physics", Elsevier Science Publishers 1986, D. Dewit 
and J. Smith. 

M. Kirson, Ann. Phys. 66 (1971) 624. 

G.E. Brown and M. Bolsterh, Phys. Rev. Lett. 3 (1959) 472. 

O. Kaczmarek, F. Karsch, P. Petreczky and F. Zantow, Phys. Lett. B 543 (2002) 
41. 

G.E. Brown, A.D. Jackson, H.A. Bethe and P.M. Pizzochero, Nucl. Phys. A560 
(1993) 1035. 

F. Karsch, K. Redhch, and A. Tawfik, Phys. Lett. B571 (2003) 67. 

P.B. Braun-Munzinger, K. Redlich, and J. Stachel, Invited review for "Quark 
Gluon Plasma 3", eds. R.C. Hwa and Xin-Nian Wang, World Scientific; 
|nucl- th/0304013( 

P. Fachini, for the STAR Cohaboration, Nucl. Phys. A715 (2003) 462. 
J. Adams et al., Phys. Rev. Lett. 92 (2004) 092301. 
R. Rapp, Nucl. Phys. A725 (2003) 254. 

W. Broniowski, W. Florkowski and B. Hiller, Phys. Rev. C68 (2003) 034911. 

G. E. Brown, J. Stachel and G.M. Welke, Phys. Lett. B 253 (1991) 19. 
G.E. Brown and M. Rho, Phys. Rev. Lett. 66 (1991) 720. 



G.E. Brown and M. Rho, |nucl-th/0305088, Phys. Repts., to be published. 



K. Langfeld, H. Reinhardt and M. Rho, Nucl. Phys. A622 (1997)620. 
V. Koch and G.E. Brown, Nucl. Phys. A560 (1993) 345. 
F. Karsch, Nucl. Phys. A698 (2002) 199c. 

O. Kaczmarek, S. Ejiri, F. Karsch, E. Laermann and F. Zantow, 
|hep-lat/0312015t Talk at Finite Density QCD, Nara, Japan, 10-12 July, 2003. 

Proc. of A RIKEN BNL Research Center Workshop, "Lattice QCD at finite 
temperature and density", Feb. 8-12, 2004, Brookhaven National Laboratory, 
USA. 



D.E. Miller, |hep-p h/000803T; Phys. Repts. in preparation. 

S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, |hep-lat/03r2C)37, 

T. Yamazaki et al. (CP-PACS Collaboration), Phys. Rev. D65 (2002) 014501. 

L.P. Kadanoff and G. Baym, "Quantum Statistical Mechanics", W.A. 
Benjamin, New York 1962. 



40 



[38] E. Shuryak and 1. Zahed, Phys. Rev. D69 (2004) 046005. 

[39] H. Stocker, private communication, to be published. 

[40] G.E. Brown and M. Rho, Phys. Rept. 269 (1996) 333. 

[41] M. Harada, Y. Kim and M. Rho, Phys. Rev. D66 (2002) 016603. 

[42] E.V. Shuryak and G.E. Brown, Nucl. Phys. A717 (2003) 322. 

[43] G.E. Brown and M. Rho, Phys. Rept. 396 (2004) 1. 

[44] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25 (2000) 1. 



41 



